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Abstract: Conventional control charts are one of the most important techniques in statistical process 
control which are used to assess the performance of processes to see whether they are in- or 
out-of-control. As traditional control charts deal with crisp data, they are not suitable to study 
unclear, vague, and fuzzy data. In many real-world applications, however, the data to be used 
in a control charting method are not crisp since they are approximated due to environmental 
uncertainties and systematic ambiguities involved in the systems under investigation. In these 
situations, fuzzy numbers and linguistic variables are used to grab such uncertainties. That is why 
the use of a fuzzy control chart, in which fuzzy data are used, is justified. As an exponentially 
weighted moving average (EWMA) scheme is usually used to detect small shifts, in this paper a 
fuzzy EWMA (F-EWMA) control chart is proposed to detect small shifts in the process mean when 
fuzzy data are available. The application of the newly developed fuzzy control chart is illustrated 
using real-life data. 


Keywords: crisp data; fuzzy data; fuzzy control charts; EWMA charts; fuzzy EWMA charts 





1. Introduction and Literature Review 


Process performance and proficiency can be enhanced by decreasing variability. This can be 
attained with the assistance of statistical process control (SPC) methods. The aim of a statistical 
process control method is to achieve stability and to improve the efficiency of a process by reducing 
the variability involved. Conventional control charts are one of the tools to achieve this objective. 
Conventional control charts, also known as Shewhart control charts, are used to detect shifts in a 
process. The mostly used Shewhart control charts when a variable quality characteristic is monitored 
include “x-bar and range” and “x-bar and s(standard deviation)” control charts. However, these charts 
are not able to detect small shifts. Another problem is that these charts are not independent of each 
other. In order to detect small shifts in the process mean or in the process variation, an exponentially 
weighted moving average (EWMA) scheme is used. Nonetheless, the construction method and the 
interpretation involved in a traditional Shewhart chart are different from the ones in an EMWA control 
chart. In other words, the plot of a Shewhart control chart shows independent data points and follows 
the law of the sampling distribution of the statistic. However, in an EWMA or in a cumulative sum 
(CUSUM) scheme, a statistic corresponding to a parameter is estimated directly. When a process 
monitored using an EWMA or a CUSUM control chart goes out-of-control, the data moves gradually to 
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the next level. Otherwise, the chart will show a discrepancy about the centerline with small variations 
when the process is in-control [1]. In addition, the EWMA and the CUSUM control charts are usually 
used to detect small shifts in a process. Yang et al. [2] proposed a new EWMA control chart to monitor 
the mean of a process with a variable quality characteristic. It can be applied in a situation when the 
distribution of the available process data is either unidentified or non-normal. Although their chart 
was shown to have a very good performance, the Arcsine-EWMA control chart developed also by 
Yang et al. [2] performed better. To monitor processes with attribute characteristics, however, the p 
and the c charts are the most widely employed schemes. 

A control chart consists of three lines; the center (average) line, the upper, and the lower control 
limits (lines). The strength of any chart can be judged by its ability to identify the shifts and the 
irregularity of a process. This strength is directly related to the choice of its control limits when the 
chart is designed. This also applies to fuzzy control charts introduced by Raz & Wang [3] for the first 
time, when they suggested the use of two approaches: a probabilistic approach and a membership 
approach. Ever since, this type of chart has gained much interest in the literature. The control charts 
are recommended by Kanagawa [4] for linguistic terms obtained by experts’ opinion and regarded 
as fuzzy data are different from the one proposed by Wang [5]. It monitors not only the process 
average but also the process variability. In order to decrease the emergence of false alarms, as well as to 
enhance the speed to identify actual errors, El-Shal [6] used fuzzy logic to amend SPC rules. Rowlands 
& Wang [7] proposed a fuzzy SPC evaluation and control method which combined the traditional 
statistical process control methodology with an intelligent system approach. Gülbay & Ruan [8] 
recommended «-cut control charts to monitor attribute data to adjust the stiffness of the inspection 
using triangular fuzzy numbers. Cheng [9] developed a fuzzy control chart to deal with a process that 
gives fuzzy outcomes. He designed two fuzzy control charts to directly observe the fuzzy results in 
order to determine the process quality. A direct fuzzy approach was introduced by Gülbay [10] as a 
substitute to fuzzy control charts. This approach can be applied on fuzzy attribute data obtained by 
vague events. To control process mean, Faraz [11] developed a fuzzy control chart with a warning line 
alongside an upper control limit. Erginel [12] presented the theoretical structure of fuzzy individual 
and moving range control charts with «-cuts using a-level fuzzy median transformation techniques. 
Sentiirk [13] designed fuzzy (X — R ) and (X — S) charts with a-cuts using «-level fuzzy midrange 
transformation techniques. They used triangular fuzzy membership functions with fuzzy numbers 
(a, b, c). Later, Sentiirk [14] proposed a fuzzy regression control chart. Fuzzy variable control charts 
deal with fuzzy variable data, for example, fuzzy measurement (length, width, height). On the other 
hand, fuzzy attribute control charts study qualitative data which cannot be expressed in numbers, 
for example, fuzzy attributes (defective, nonconformities). While some fuzzy variable control charts 
include (X — R ) and (X — 5), fuzzy attribute control charts include (j,¢, i). If the quality-related 
characteristics cannot be represented in arithmetic form, such as characteristics for appearance, softness, 
color, and so forth, then control charts for attributes are used. Product units are either classified as 
conforming or nonconforming, depending upon whether or not they meet specifications, or the number 
of nonconformities (deviations from specifications) per unit is counted. There are four major types 
of control charts. The first type, the x-chart (and related x-bar, r-, and s-charts), is a generic and 
simple control chart. The x-chart is designed to be used primarily with variable data, which are 
usually measurements, such as the length of an object, processing time, or the number of objects 
produced per period. The other types of control charts are the p-charts, which are used with binomial 
data, the c-charts, which are used for Poisson processes, and the u-charts, generally designed for 
counting defects per sample when the sample size varies for each inspection. According to Kaya & 
Kahraman [15] when the value is given as “around” and “approximate”, a triangular fuzzy number 
is mostly used, and when the value is given as “between”, two trapezoidal fuzzy numbers are used. 
The examples of a triangular fuzzy number and trapezoidal fuzzy numbers are approximate speed 
(20 km) and speed between (25 km to 30 km). If the vertex of a trapezoidal fuzzy number has taken a 
common value, it is transformed into a triangular fuzzy number. 
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Sentiirk et al. [16] introduced fuzzy 1 control charts to monitor fuzzy attribute data. Kaya [15] 
designed fuzzy control charts for fuzzy measurements of the related quality characteristics. A fuzzy 
multivariate exponentially weighted moving average (F-MEWMA) control chart was suggested by 
Alipour [17] to be implemented on a process in the food industry. Sentiirk [18] suggested a fuzzy 
EWMA (F-EWMA) control chart to study fuzzy univariate process data using «-cuts. The design 
of this scheme is based on using triangular fuzzy numbers and the fuzzy midrange technique to 
study fuzzy data by converting it into scalar form. Their future recommendation was to apply the 
same technique on conventional control charts to make them able to work in fuzzy environment. 
Khademi [19] developed two fuzzy approaches to design control charts observing fuzzy averages 
and ranges. These approaches that are based on fuzzy mode and fuzzy base rules are better in the 
sense they do not use a defuzzification method and are more informative than the existing fuzzy 
control charts. They demonstrated the performance of their charts with the help of a real-life example. 
Recently, Kahraman [20] introduced two fuzzy control charts to monitor fuzzy outcomes in order to 
see whether a process is in-control or out-of-control. It should be mentioned that the capability of a 
process can also be measured with the help of a fuzzy control chart to check whether a process is in- or 
out-of-control [15]. 

The recommendation provided by Sentiirk [18] is utilized in the current paper on another 
conventional scheme proposed by Yang et al. [2] to design a new fuzzy EWMA chart that can be 
used to detect small shifts on the mean of a process in a fuzzy environment. This makes the newly 
design chart more applicable, more flexible, and more informative in real-life quality control problems. 
It gives practitioners an opportunity to study intermediate values of fuzzy parameters that is not 
possible with crisp parameters in a fuzzy environment. Although Yang et al. [2] used approximate 
values in their example, they are treated as fuzzy numbers in the current work, based on which a 
fuzzy form of his chart is proposed. This newly developed fuzzy chart will be applied on real-life data 
collected from the food industry. 

The structure of the remainder of the paper is as follows. A brief background on the fuzzy 
variable and attribute control charts and the fuzzy transformation approaches is given in Section 2. 
The conventional EWMA chart proposed by Yang et al. [2] is described in Section 3 in detail. 
The fuzzy EWMA chart to monitor the process mean in a fuzzy environment is designed in Section 4. 
The application of the proposed chart is investigated in Section 5 using cooking oil filling data. Finally, 
the conclusion and recommendation for future studies are given in Section 6. 


2. The Fuzzy Attribute and Variable Schemes 


Fuzziness and ambiguity may have several reasons including lack of knowledge, chance, 
imprecision, incompetence to perform adequate measurements, and inborn fuzziness in the process. 
Fuzziness in humanistic systems is dealt with in a mathematical way with the help of fuzzy sets [21]. 
The fuzzy set theory provides an opportunity to study human subjectivity in a scientific manner with 
the assistance of fuzzy sets [22]. 


2.1. Fuzzy Transformation Approaches 


When fuzzy data are used, it is essential to represent the fuzzy sets linked with the linguistic data 
in the sample by some scalar (transformation) numbers for additional calculations. There are four fuzzy 
measures of central tendencies in the literature used to calculate transformed numbers. These measures 
include (1) fuzzy mean, (2) fuzzy median, (3) fuzzy mode, and (4) fuzzy midrange ([3,5]). 


2.2. Some Fuzzy Formulae 


The formula to compute the fuzzy average is [21]: 


free = Sao Ur (X) 
a Ur(X) 
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To determine the fuzzy median, the curve under the membership function is divided into two 
halves. As a result, the fuzzy median becomes 


Í. med 
Ímedian = f Ur(X)dx. 
The fuzzy mode is the value of the base variable where the membership function equals 1. Thus, 


fmode = {X/Ugf(X) = 1}. 


The fuzzy mid-range is the mean end points of the «-cut as 


1 
Ff mr = za + Ac). 


According to [18], there are no solid reasons to decide which technique is better to use, but most 
researchers applied the a-level fuzzy midrange technique in their works due to its simple form 
of formula. 


2.3. Triangular Fuzzy Number 


A triangular fuzzy number A is a fuzzy number with its membership function defined by three 
numbers a, < a2 < az where the base of the triangle is the interval |a,,a3] and the vertex is at X = ap. 


Example. In the following Figure 1, there is triangular number defined by three numbers (a1 = 2, az = 3, a3 = 4), 
where the base of the triangle is the interval |2, 4| and the vertex is at X = 3. 
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Figure 1. Triangular Fuzzy Number. 


The basic forms of the variable and the attribute fuzzy control charts were given in [20]. 
The attribute control charts are p-chart, np-chart, c-chart, and u-chart, while the variable charts 
include (X — R) and (X — S) schemes. Fuzzy forms of these charts make them more flexible, sensitive, 
and informative. A fuzzy variable quality characteristic in m samples, each with the size of n was 
represented using a triangular fuzzy number denoted by X; = (Xia, Xiv, Xic), where i = 1,2,3,...,m. 
In addition, each sample in a fuzzy attribute chart, for example, the np control chart developed to 
study the number of nonconforming units in a fuzzy environment, was presented by a triangular fuzzy 
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number (dia, dip, dic) with the mean of (np, npp npe), where np, = Liz fu nD, = Liz, fi 


Lir fie, Hence, the control limits of the fuzzy np chart are [20]: 


UCL = (nP, +3 \/ np, (1 E Pa) Pp +3 y np,(1 E Pp), np, + 3y/np_(1 ~~ P.)) (1) 


CL= (nPar "Py np.) (2) 


LCL = (np, —3y/np,(1 — Pa): Py — 3y/ np, (1 — Py), np, — 3y/ np. — P,)). 6) 


2? and np, = 





3. The EWMA Control Chart to Monitor Process Mean 


A new version of the conventional EWMA control chart, initially presented by Montgomery [23], 
was proposed by Yang et al. [2] to monitor small shifts in a process mean. The details of this chart are 
given as follows. A random sample of size n including the random variables X1, X2, X3, . . -, Xn is taken 
from a process with a mean y in subgroup i. Let 


Yke J=, rerh L e (4) 


In a subgroup, define a Bernoulli random variable as 


1; Y;>0 
I; = 7 ] : 5 
] l 0; Otherwise (5) 


Let S be the total number of Yj > 0. Then, S = X; Ij would follow a binomial distribution 
with the parameters (n, p) for an in-control process with p = P{Y; > 0} as the probability of having 


a positive difference (Xj; — X) > 0, in which X is the estimate of the process mean u. Although 
the S statistic in Yang et al.'s EWMA chart [2] follows a binomial distribution, the difference with 
the conventional np chart is that the binomial variable is not the count of non-conforming units in 
the sample, rather the number of X; values in a sample that are above the in-control process mean. 
This chart is designed to detect small shifts rapidly and effectively. 

The EWMaA statistics in Yang et al.’s EWMA chart [2] is defined as 


EWMAg; = AS; + (1—A)EWMAsz 0<A<1 (6) 


where S; represents the value of S in the ith subgroup. The initial value of this statistic is EWMAg; = np. 

As the mean of the statistic is E(EWMAs;) = np and the variance is Var(EWMA<g;) = np(1 — p) with 

the limiting value of Var(EWMAg;) = z [np(1 — p)], then depending on the number of the subgroup 

in which a sample is drawn, the EWMA chart proposed by Yang et al. [2] is designed by the control 

limits in Equations (7)-(9) and (10)-(12) for a large and a small number of subgroup, respectively. 
When the subgroup number is large: 





UCLewmas = np hy] noc -= p)] (4) (7) 


CLEWMAs = p (8) 





LCLewMaA, = "p — efra — p)] (545) . (9) 
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When the subgroup number is small: 














UCLewma, = 1p + kylnva p)] (z a i) [1— (1- A)*] (10) 
CLEwMA, = "p (11) 
LCLewmMas = ^p — njma p)] (; 1 z) [1—(1—A)*]. (12) 


This version of the EWMA chart is used in the current work to design a fuzzy EWMA chart in 
Section 4. However, before proposing the chart, a brief background is provided in the next section to 
describe an existing fuzzy EWMA chart. 


The Existing Fuzzy EMWA Control Chart 


In the fuzzy EWMA scheme proposed by Sentiirk [18], z; is the tth exponentially weighted 
moving average, X; indicates the sample mean, 0 < A < 1 is a real constant, X is the overall sample 
mean, m is the number of samples, t = 1,2,3...m is the sample number, and zo = X. Assuming Xj; 
to be independent random variables with a known variance Lae the variance of z; becomes 07; = 
2 (4,)-( A)”], where n. is the sample size. Then, when t is small the control limits in the 
traditional EWMA control chart are as follows: 

















- 2 
CLewMa = X (14) 
= 2 
LCLewma = X | (2 SE (1— ay]. (15) 


However, as t increases o7,; tends to a limiting value Ta = a ( z). As such, the control limits 
in the traditional EWMA control chart become 


= o2 À 
Wei = es “(44) (16) 
CLewma = X (17) 
LCL zaia i ae (18) 
EWMA = Fike 


On the other hand, the control limits of the fuzzy EWMA scheme proposed by [18] when ø is 
known and t is small are 





UCLEewMa = Coe rA ag T Cno, cal (z5) 1- Oe a)” (19) 
ČLewma = (Kauko x) (20) 





UCLewMa = (Xa, Xb, Xe) — Sp Cuon aJ) [1 — (1-14)”]. (21) 
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However, when t is large, the control limits become 








~ = = = 3 À 
UČLewma = (Xa Xy Xe) + T Onon] (75) (22) 
CLewma = (Xa, Xp, Xc) (23) 
UCL san Je ae) de (24) 
EWMA — ar br Ac Vn arUb, Uc 2—1)’ 
In case g is unknown and t is small, the control limits are 

~ = = = Bt Ed: as À 
UCLewma = (Xa Xp, Xc) + RaRo RA (z2) 1- (1-14) (25) 
CLewma = (Xa, Xo, Xo) (26) 

x Sa is EE À 
UCLewMa = (Xa, Xp, Xe) T 3(Ra, Ry, rJ) [1 = (1 = à)”] . (27) 


Finally, when g is unknown and t is large the control limits are 


~ = i: = a À 
UCLewmMa = (Xa, Xp, Xc) + 3(Ra, Ro, Re) (545) (28) 
CLewma = (Xa, Xb, Xc) (29) 

x = = 5 =c À 
UCLewma = (Xa, Xp, Xc) 7 3(Ra, RBI (74) . (30) 


Note in Equations (25)-(30) that Ra, Rp, and R, are ranges. 


4. The Proposed FEWMA Scheme Based on the np chart 


Based on the background provided in Section 3, the control limits of the new fuzzy EWMA 
scheme to detect small shifts in the process mean when fuzzy process data are available are proposed 
as follows. Similar to the limits shown in the existing fuzzy EWMA chart, assuming an unknown 
standard deviation of the process, the case that mostly happens in practice, the limits are provided 
depending on whether t is small or large. 


4.1. æ Is Unknown and t Is Small 


The control limits, in this case, are proposed as 











UCLewanas = ("PaPa 1P) +Ky/ (32x) 0 (1-A) apC- Pa) PO- Pa) Pe -PA (89) 


CLEWMAs = ("Py "Py, npe) (32) 











LCLewnas = (nPar Pp, NPs) Ky (4s) 0 (1—A)"'][mp, (1 — Pa)l, [Py (1 — By)], nP. — ,)] (33) 
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which can be written in the following forms: 








UCLewmas = (mPa + Ky) (oA) f= (1 -APPC Pade MP 


FK| (aèr) E- 1 APC = Pa) oP, 














CLEWMAs = (np, "Py, NP.) 








iol er (up, Ky/(s4z) 0 (1—A)*][np, (1 — Pa)l, "By 
Ky (2x) E- (1 = A)", = Pade mB 
K (4x) I - 0- 0)*9i0p. 4-7.) 








8 of 13 


(34) 


(35) 


(36) 


In Equations (34) and (36), the parameter K is chosen such that the chart exhibits a desired 


in-control average run length. 


4.2. o Is Unknown and t Is Large 


The control limits in this case are 








Het Hae Punpa +K | (52) P.O -POL PC - Pol P.A -PAN 


CLewMag = ("Pq "Pp, nP) 








ECL swans = (oy nb P)-Ky] (<2, ) rp. — Pa) PoC- Pad med — Pe) 


which can be rewritten as 











CLewmas = ("Pq "Pp MP) 
LCLewaas = (Pa -KY GAD P,O — Pad] MP 
- K/A) inp, (1 = Py], np. 
-Ky (32x) [np -(1 = B)] )- 











(37) 


(38) 


(39) 


(40) 


(41) 


(42) 


Then, depending on the fuzzy measures of central tendencies used for the fuzzy transformation 


approach discussed in Section 2.2, the following limits are obtained. 


4.3. The Control Limits in the -Cuts Fuzzy EWMA Control Chart 


The control limits in this fuzzy transformation approach are 
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UCT*ewaas = (HP, +Ky z2 (np — )), nP," 























+K y 4s (0P, (1 — B,*), mB, (43) 
+K y (nP (1 - ay 
CL*EWMAs = (np,", np," np") (44) 
LCL" tewma; = (np," +Ky/ 45 ey CT PC — Pa“) np," 
HKP — By") mB" (45) 





+Ky/ (np — np,"))). 


4.4. The Control Limits in the «-Cut Fuzzy Median EWMA Scheme 


The control limits of the proposed fuzzy EWMA scheme when this type of fuzzy transformation 
approach is used are shown in Equations (46)—(48). 














UCTS ed —EWMAs = $ (Pa + MP py" + MB.) +4K y/ xy (mp ,"(1 — By*)) + OPC — By") + (mp1 — mp) (46) 





ee: ma 7 7 
CL* med—EWMAs = 3 ("Pa + np," +np,") (47) 

















ECL nedewM Aas = § (Py + By + nP) -IK 94x (MP8 Pat) + P — By") + (BAL — mp") (48) 


5. The Application 


In this section, the application of the proposed fuzzy EWMA chart is demonstrated using 
fuzzy observations obtained from a cooking oil filling process involved in the food industry in 
Pakistan. In this example, fuzzy data are collected from the process, where fuzzy observations are 
treated as triangular fuzzy numbers (Xj, Xpj, Xci) in each sample taken at different times. In addition, 
the proportion of the values above a given threshold (zero) will be treated as a triangular fuzzy number 
(Pair Pois Pci). The data consist of 10 samples, each with 10 fuzzy observations, and the proportion is 
obtained from each sample based on the given threshold. 

In this application, the chart parameters used are K = 2.58, A = 0.2, n = 10, pa = 0.35, pp = 0.48, 
« = 0.65, and pe = 0.64. Here, the parameter setting K = 2.58, A = 0.2, and P = 0.65 is obtained using a 
simulation approach described in the next section. In addition, pa = 0.35 shows the proportion of “a” 








above its given mean, p, = 0.48 implies the proportion above the mean of “b”, and p; = 0.64 indicates 
the proportion above the mean of “c”. This is done due to the fact that each value is being used as a 
triangular fuzzy number. The control limits are calculated as 


UCLewma, == (4-75, 6.11 7.66) 


CLewma, == (3.5, 4.8, 6.4) 
LCLewma, == (2.22,4.46,5.13). 


Then, the limits of the a-cuts fuzzy EWMA control chart are obtained based on Equations (43)-(45) as 
UCL* waa, == (5.69, 6.15, 6.7) 


LCL* ew, == (2.99, 3.44, 4.0). 


However, when the «-cut fuzzy median EWMA scheme is used, the control limits based on 
Equations (46)-(48) are obtained as 
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——_ 


UCL” med—EWMAs = 5.62 
CE" med -EWMAg = 4.84 
LOE ned EWMAg = 4.05. 


When the «-cut fuzzy median EWMA chart is employed on the first sample we have 





MPa at = nNPa1 F a(nppı = NPa) (49) 
Np" ci = npea + a(npyr — npea) (50) 
1 : 
5  med—EWMAsj = 3 ("P"aj + Pej + MP%cj)i j= 1,2)... - (51) 


Then, at « = 0.65 we have 
np 6 = 4.65 


As a result, the statistic in the first sample becomes 


a Sedat: 
S* med-EWMAs1 == 9 


Now, since the statistic falls within the control limit, i.e., 4.05 < 5 < 5.61, the process is in-control 
based on the first sample drawn. 

The proportions of the fuzzy observations above the process mean are first calculated for each 
sample in Table 1. The detailed calculation using the R statistical software is given in Appendix A. 
Then, the control limits of the proposed fuzzy exponentially weighted moving average control chart 
using the a-cut median approach are computed. The results show that the process is in-control. 


Table 1. Monitoring the fuzzy oil packaging process using the proposed fuzzy exponentially weighted 
moving average (EWMA) scheme. 


Sample Fuzzy Proportion a-Cuts for the Proportion Fuzzy EWMAs Fuzzy (EWMA-med) 4.054<S* pmwas"<5.62 





1 4,5,6) (4.65,5,5.35) (3.6,4.84,6.32) 4.92 In-contro 
2 3,4,5) (3.65,4,4.35) (3.8,4.8,5.8) 4.80 In-contro 
3 5,6,7) (5.65,6,6.35) (3.4,4.4,5.4) 4.40 In-contro. 
4 2,3,4) (2.65,3,3.35) (5.05,5.4,5.75) 5.40 In-contro 
5 6,5,7) (6.65,6,6.35) (4.80,6.3,6.46) 4.60 In control 
6 3,4,5) (3.65,4,4.35) (4.6,5.6,6.6) 4.93 In-contro 
7 4,6,7) (4.65,3,2,3.25) (3.40,4.6,4.7) 4.63 In-contro 
8 3,4,5) (3.65,4,4.35) (4.6,5.6,6.6) 4.93 In-contro 
9 4,5,6) (4.65,5,5.35) (3.2,4.2,5.2) 4.20 In-contro 
10 6,7,8) (6.65,7,7.35) (4.4,5.4,6.4) 5.40 In-contro 








The Simulation 


In this section we simulated the data using MATLAB code to calculate Fuzzy average run length 
(ARLI1s) when the process is not in-control shown in Table 2. Similarly, conventional ARL1s are 
presented in Table 3. Comparison shows that the fuzzy ARLIs are less than the conventional ARLIs. 
This means that the fuzzy control scheme is quicker in indicating small shifts in the process as compared 
to the conventional chart. This is the strength of the proposed fuzzy control chart compared to the 
conventional chart. 

The following algorithm was used to find the control chart coefficients and the ARLs. 


1. Generate a binomial random variable. 
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2. Determine the in-control limits using K = 2.58, P = 0.65, A = 0.2 so as the in-control ARL, ARLo, 
as a function of K, P1, and A, becomes 371 using 10,000 simulations, each with 25 samples. 

3. Similarly, to estimate the ARL1, the process is simulated 100,000 times for an out-of-control 
process with P1 = 0.65. In this case, the ARL1, which is a function of K, P1, and A, becomes 2.34. 

4. The same simulation process is repeated in a fuzzy environment, the difference being the use of 
fuzzy random variables. 


Table 2. Fuzzy average run length (ARL1s) of the proposed chart when A = 0.2, K = 2.84, P = 0.65. 








P 0.25 0.35 0.45 0.55 0.613 0.65 0.75 0.85 0.95 
n 

9 1.0 1.0 1.0 1.00 3.61 5.57 3.04 1.02 1.0 
10 1.0 1.0 1.0 1.55 4.11 3.06 3.45 1.06 1.0 
11 1.0 1.0 1.0 4.70 3.61 3.48 3.12 1.05 1.0 
12 1.0 1.0 1.0 1.37 1.37 3.15 3.06 1.30 1.0 
13 1.0 1.0 1.0 1.02 1.50 2.82 2.37 1.23 1.0 
14 1.0 1.0 1.0 1.005 1.13 2.99 2.87 1.05 1.0 
15 1.0 1 1.003 3.07 4.26 16.91 2.26 1.01 1.0 


Table 3. The ARL1 of the existing chart when A = 0.2, K = 2.84, P = 0.65. 








P 0.25 0.35 0.45 0.55 0.613 0.65 0.75 0.85 0.95 
n 

9 1.0 1.0 1.34 5.41 48.06 188.53 6.35 1.32 1.0 

10 1.0 1.0 1.14 5.40 86.40 330.32 6.17 1.07 1.0 

11 1.0 1.0 1.23 4.26 105.68 575.78 5.48 1.04 1.0 

12 1.0 1.0 1.11 4.54 101.57 743.43 4.56 1.08 1.0 

13 1.0 1.0 1.21 4.41 104.52 1000.35 4.64 1.10 1.0 

14 1.0 1.0 1.05 3.32 181.62 1200.14 3.79 1.0 1.0 

15 1.0 3.33 1.03 3.21 178.23 1330.43 3.87 1.0 1.0 


6. Conclusions 


Statistical process control (SPC) techniques play a very important role to effectively utilize 
manufacturing resources and to enhance the quality of the products manufactured. Statistical quality 
control charts are one of the most important techniques in SPC. Among these schemes, the exponentially 
weighted moving average (EWMA) control charts have been shown effective in detecting small process 
shifts. While precise data are required in the available conventional control charts, fuzzy data are 
generated due to uncertainties involved in the practice. Thus, the use of fuzzy control charts are 
justified, as they are more sensitive and have more flexibility compared to the conventional charts. 
In this paper, a novel fuzzy EWMA control chart was proposed to monitor processes with fuzzy data. 
The newly developed FEMWA control chart was applied on fuzzy data collected from the cooking oil 
industry and showed effectiveness. The same analysis can be extended in the future when, instead of 
rectangular, trapezoidal fuzzy numbers are used. In addition, a comparative investigation is required 
between the performances of the conventional and the fuzzy EWMA schemes in terms of the average 
run lengths (ARL) criterion. 
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Appendix A 


Calculation with R 

rm<(list=1s()) 

n=10 

Pa=0.35 

Pb=0.48 

Pc=0.64 

K=2.58 

#Alpha-cut fuzzy EMWA control chart (OK) 

#Var=n*p*(1-p)*0.333 #where a/2-a=0.3333,a=0.2 

Pacut=Pa+0.65*(Pb-Pa) 

Pb=Pb 

Pccut=Pc+0.65*(Pb-Pc) 

CLMedian=(n*Pacut+n*Pb+n*Pccut) /3 
UCLmedian=CLMedian+(1/3)*2.58*0.33333*sqrt(n*Pacut*(1-Pacut)+n*Pb*(1-Pb)+n*Pc*(1-Pc)) 
LCLmedian=CLMedian-(1/3)*2.58*0.33333*sqrt(n*Pacut* (1-Pacut)+n*Pb*(1-Pb)+n*Pc*(1-Pc)) 
data.frame(LCLmedian,CLMedian,UCLmedian) 








controlimits=c(LCLa=n*Pa,CLb=n*Pb,UCLc=n*Pc) 








X1=c(4*0.2+0.873.5,5*0.2+0.874.8,6*0.2+0.876.4) 
mean(X1) 
X2=c(3*0.2+4*.8,4*.2+5*.8,5*.2+6*.8) 
mean(X2) 
X3=c(5*0.2+3*.8,6*0.2+4*.8,7*.2+5*.8) 
mean(X3) 
X4=c(2*0.2+5*.8,3*0.2+6*.8,4*.2+7*.8) 
mean(X4) 
X5=c(5*0.2+2*.8,6*0.24+3*.8,7*.2+4*.8) 
mean(X5) 
X6=c(5*0.2+2*.8,6*0.2+3*.8,7*.2+4*.8) 
mean(X6) 
X7=c(1*0.2+5*.8,2*0.2+6*.8,3*.2+7*.8) 
mean(X7) 
X8=c(3*0.2+1*.8,4*0.2+2*.8,5*.2+3*.8) 
mean(X8) 
X9=c(4*0.243*.8,5*0.2+4*.8,6*.2+5*.8) 
mean(X9) 
X10=c(6*0.2+4*.8,7*0.2+5*.8,8*.2+6*.8) 
mean(X10) 
data.frame(mean(X1),mean(X2),mean(X3),mean(X4),mean(X5),mean(X6),mean(X7),mean(X8), 
mean(X9),mean(X10)) 
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